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MAGNETOPLASMA DIFFUSION AT F2 REGION ALTITUDES 


P. C. Kendall 


ABSTRACT 

The magnetoplasma diffusion equation of the F2 region of the 
ionosphere is derived in its primitive form, allowing for electric 
currents flowing between conjugate points of the dynamo region, and 
for temperature variations. 

It appears that the standard view of the ambipolar diffusion 
problem is correct. The ion-electron pairs may be thought of as 
sliding down the magnetic field lines with diffusive motion, while 
the whole field line drifts at right angles to itself. 

The electric current along a field line would affect the ioniza- 
tion only if its magnitude were 100 times greater than Dougherty's 
estimate 1 of 10 ~ 7 amps/m 2 . 
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MAGNETOPLASMA DIFFUSION AT F2 REGION ALTITUDES 


P. C. Kendall 


INTRODUCTION 

In this note the magnetoplasma diffusion equations of the F2 layer are de- 
rived in an attempt to lay down a solid foundation for future work. This is nec- 
essary for two reasons. First, it is clear that many future computations will 
be based upon current theory, which should therefore be refined as far as pos- 
sible. Second, there has been much recent discussion, as evidenced by the work 
of Chandra 2 , Kendall 3 , Ferraro 4 , Varnasavang 5 and Gliddon 6 . 

The possible effects of a hypothetical electric current nearly along the mag- 
netic field lines will be included because Dougherty 1 has suggested that such an 
electric current may arise from asymmetry in the dynamo region. 

The electron and ion temperatures, T e and Tj, will be retained as different 
functions of position throughout, so that even in the case when there is no elec- 
tric current along the field lines, the equations developed will be as useful as 
possible. The acceleration of the plasma will be ignored. 

EQUATIONS OF THE PROBLEM 

Suffixes i , e and n will be used to distinguish quantities associated with 
the ions, electrons and neutral atoms. For convenience, it will be assumed that 
at F2 region levels (300 km) the neutral atmosphere is 0 and the ions are 0 + . 

The plasma will be assumed to be electrically neutral, so that if N denotes the 
electron density, 


N = N. = N e . (1) 

Denote by m the mass, by v the velocity and by p the partial pressure of a par- 
ticular constituent; by g the acceleration due to gravity with 

|g| = 890 cm/sec 2 , (2) 


by e the electronic charge 

e = 1.602 x 10" 20 e.m.u. 


( 3 ) 
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by B the magnetic induction vector 

I B | ~ 0.3 gauss , (4) 

by £ the density of electric current 

l = eN fe ~Ze) ( 5 ) 

and by v. , v . e , v ia the collision frequencies. At F2 region levels 

v in ~ 1 sec -1 , v en ~ 35 sec -1 , v gi ~ 1500 sec -1 . (Q) 

Taking m n = m. and assuming that m e << m. , the equation of motion of the ion- 
electron gas is 

Nm e v en (v. -v n ) + \ Nm. v in (vj -v n ) 

= ' V(p e + Pi ) + Nm. g + J x B (7) 

Using Equation (5), Equation (7) may be transformed to either of the two identical 
forms 

T N ( m i ^in +2m e V en) ( $ ~ Zn) 

" “ ^P e + Pi ) + Nm i J| + ~ l + l X 2 (8) 

and 

T N ( m i v in + 2m e v en ) (v e -v„) 

m i V in 

= - V(p e + Pi ) + Nm. g + 2e“ J + I x 5 * (9) 
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Denote by a 11 the component of an arbitrary vector a along a field line in the 
same direction as the magnetic field and by a the component of a at right angles 
to a field line such that a 1 = a - a 11 . For future use the component of Equation 
(8) along a field line will be required along with the component of Equation (9) at 
right angles to a field line. Thus 

T N K y in + 2m e ^en) ($" " &") 

rn V 

= " ? (p e + Pi ) + Nm i Ji" + “ J" ( 10 > 


and 



bn 


+ 2m e 



i . i m i bn 1 1 

= " V P e + Pi ) + Nm i 8 — J + J X B . 


(li) 


Equations (10) and (11) together constitute the equation of motion of the ion- 
electron gas . The equation of motion of the electrons is 


Nin ( v “ v 1 + Nrn v . ( v ~ v. ) 

£ Gtl \ / 0 G 1 y ^ / 


- - V p - Ne(E + v x B ) , 

r\j C \ /X, ryj* '-Xj / 


( 12 ) 


where E denotes the electrostatic field. 

Denote by E t the electrostatic field related to motion of the electrons, by 
E 2 the electrostatic field due to pressure gradients and by E 3 the external elec- 
trostatic field. As the geomagnetic field is supposed to be unchanging, these 
electric fields are supposed without loss of generality to arise'from potentials 
Dj , fl 2 and ft 3 , where 


E 


l 


- Vfi 


i ’ 


vn. 


- vn. 


(13) 
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and if s denotes the arclength and T denotes the unit tangent along a field line 


E 


1 


III 

e e n 




+ 


rn 


e 




( 14 ) 


and 


E 


2 



T Bp e 

eN 3s 


giving 


(15) 


n 


1 




and 


n 


2 



J_^Pe 
eN 3 S 


ds , 


(16) 


(17) 


where v 11 = v 11 T for all constituents, and the level s = 0 is conveniently chosen. 
Then flj and 0 2 are known in terms of the other variables, and Equation (12) 
shows that 


30 


os 


= 0 


(18) 


This shows that as defined here, the external potential 0 3 is constant along a 
field line. Note that functions F and G, arbitrary apart from satisfying the 
equations 


3F 3G 

3s ~ °- 3s " °> (19) 


may be added to O x and 0 2 . However, it is clear that these could be absorbed 
in 0 3 . It follows that the definitions of internal and external electrostatic fields 
are a matter of convenience only. From Equation (12), all three potentials 
together cause a velocity v 1 of the electrons given by 
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( 20 ) 


B + 


e e n 


= - 


eN 


7 A P. - 


m v. m v 

_1 ii T 1 + _£ 2IL i 

q 2 e r\j n 


The component of Equation (12) along a field line (i.e., Ohm's Law) has already 
been used in Equation (14). Taking m e = 9. 108 x 1CT 28 gm , Equations (3) and 
(6) show that 

m e ^en/eB « 1 • (21) 

Thus, to a high degree of approximation, 


m e ^ei 


NB 2 


T X x B 


nr v M 

e e n 


eB 2 




- B 


(22) 


In fact, the last two terms in this equation are also likely to be small. However, 
it is preferable at this stage not to jeopardize the argument by using approxima- 
tions which are not immediately justifiable. 

The system of equations has now been reduced to two Equations (11) and (22) 
in the vector quantities and v e x . 

ANALYSIS OF THE EQUATIONS 


Substituting back into Equation (11) from Equation (22) gives an equation 
for J 1 in the form 


m. v. 

i in 


m e ^ei 


6 - + i 1 " 2e 2 B 2 ( m i Vin +2me * 5 


N 


" T ( m x "in + 2m e ^en) 


B I 1 


vV- 


m v 

e e n I 

7"B X v n 


+ ^ (Pe + Pi) ‘ Nm i % 


(23) 
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Taking nij - 13. 38 x 10 24 gm, Equations (3) and (6) show that 


nr 

e e n 


eB 


m i V in 

<< — eT « 1 


(24) 


Using this result in Equation (23) it is seen that all terms on the right hand side 
are negligible in comparison with £ x B. Thus, to a high degree of approxima- 
tion, vector multiplication by B using the expansion formula for triple vector 

'Xj 

products, gives 


Nm. 


I 


V in 


2B 2 


in z/ 
e en 


v n - V Q + 


Ne 


V' 


Pe " jg X V„ 


+ 






(25) 


The corresponding relative velocity of ions and electrons at right angles to a 
field line is obtained from Equation (5) and is 




m i "in 


2eB 2 


m v 

e e n 


11 11 
v - V Q + Tj” V p 


+ 




(26) 


1 

The term V fi leads to electrodynamic drifts which are believed at most to be 
of the order of 20m/sec (Maeda 7 ). The neutral air speed | v | is unlikely to 

'%>ri 

exceed this value. Moreover, the pressure gradients will each be at most of the 
same order of magnitude as the gravitational force. It follows that as far as 
orders of magnitude are concerned, the relative velocity at right angles to a 
field caused by the electric field is 


, , , 1000m i n 

“ £ I eB " 2 ‘ 8 cm/sec , 


(27) 


while the relative velocity caused by the mechanical forces is 

m i 6 

l^e 1 " Zi I ~ "eB - ~ 2. 5cm/sec. 


(28) 
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It is clear that these relative velocities are both small compared with the esti- 
mated drift speed of 20m/sec and the estimated gravitational terminal speed 
g^* 1 ~ lOm/sec of a falling ion-electron pair. Therefore, if the electron drift 
velocity is even as high as 1 m/ s the ion velocity at right angles to a field line 
may be taken to be approximately the same as v e . Neglecting terms which again 
prove to be small, Equation (22) gives 






(29) 


This well known result shows that plasma in the F2 layer can move at right angles 
to the magnetic field lines only when an electric field is present at right angles to 
the field lines. Thus, so far, Equations (14), (18) and (29) have shown that the 
following electric fields and simultaneous drift motions may be present in the 
F2 region. 

An Electric Field Arising From or Causing Electric 

Currents Along the Field Lines 

This electric field Ej = - Vf^ is given by Equation (14). If there is an 
imposed electrostatic field of this nature from the dynamo region it will be caused 
by a difference in potential between the Northern and Southern conjugate points. 

It will appear later that the effects of ^f 11 caused by E 11 are only apparent in 
Equation (10). As m e v en << m. v. this implies that the only currents of this 
nature which might affect the electron density are such that 



or 



(30) 


so that to be consistent with div £ - 0, 


div(NvJ') ~ 0 . 


(31) 


Denote the equation of a magnetic field line by 


a(r, 6 ) ~ constant , (32) 

where r , 6, 4> are spherical polar coordinates whose axis is the axis of symmetry 
of the magnetic field. Then if B(r, 6 ) denotes the magnetic field strength, since 
div B = 0, Equation (31) gives 
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¥(a, 4>) B(r,6>)/N 


(33) 


where *P is an arbitrary function of a and 4> which is constant along the magnetic 
field lines. At this stage, it is convenient to use the direct conductivity cr along 
the field lines, giving 


E* 1 - - ecr 1 W( a , 4>) B( r , 6 ) 

r\j 1 


(34) 


Then 


1 f S 

V [cr” 1 a , cfi) B( r , a ) ] ds . 

^ Jo 


(35) 


The corresponding drift is then 


E- x B 


B 2 


x ^ J [eo-" 1 f( a, 0) B( r , 6) ] ds 


(36) 


As the primary cause of the currents lies in the dynamo region, the magnitude of 
the current depends crucially on the processes occurring there. Dougherty 1 has 
estimated that the current density along a magnetic field line may be of the order 
of 10 -7 amp/m 2 , which is 10" 12 e.m.u. Thus, so far as orders of magnitude are 
concerned 


e'l'B ~ 10 -12 e.m.u. 


(37) 


It is also clear that the function 'W3 will vary only slowly with position, whereas 
the conductivity o- is likely to vary more rapidly. It follows that 


B 


B 2 


eW 


B x V 


T 

•'O 


ds 


The order of magnitude of this can be immediately found in the case wheno- - O’(z) 
is a function of only the height z . If n denotes a local coordinate at right angles 
to a field line, such that, locally, z = s sin I + n cos I we have 


x B 
B 2 


eW 
B 2 


V 


f 


dz 

cr sin I 
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where z e denotes the height of a field line at the dip equator. Thus 


E. x B 

B 2 


eW 


cot I _ cosec I 
CT ( Z ) " CT ( z e ) 


For functions c( z) which increase with height fairly rapidly, this shows that the 
above drift velocity for a particular field line is zero at the dip equator and in- 
creases away from the dip equator. At moderate latitudes 


E. x B 

B 2 


e'P 10' 12 

~~~ "oB~ cm/sec • 


(38) 


It seems likely that this estimate of the drift due to electric currents along the 
field lines will apply to most functional forms of o -. Also, in Equation (38) it 
should be noted that cr is the local value of the direct conductivity, giving the 
local drift of this kind. 


Taking B ~ 0.3 gauss and a ~ 10' 8 or 10" 9 e.m.u. Equation (38) shows that 
the electrodynamic drift arising from E x is negligible, even for currents 100 
times larger than those predicted by Dougherty 1 . 

However, there is still a component of electric current J 11 in Equation (10), 
and this will be discussed later. 

A Pressure Gradient Electric Field E 2 = ~^ 2 

This electric field is given by Equation (15), but the combined drift caused 
by both electric field and pressure gradients is given by Equation (29) to be 


{h + eN Z p e) 


B 

B 2 


Using Equation (17) and assuming that the operators V 


1 and f 

Jo 


(39) 


ds commute gives 


i 1 i 
h + ~eN"V p e 


r k i 

% 

l 


eN ' 
J o 


-TT’ N j 


ds 


(40) 


where k is the Boltzmann constant, 


k = 1.38 x 10' 16 erg/deg . 


(41) 
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It will be sufficient to ignore the curvature of the field lines in estimating the 
magnitude of this expression. Thus, if n denotes the meridional unit normal to 
a field line, denotes the azimuthal unit vector and I denotes the dip angle 


and 


giving 


3 3 cos I 3 

3s - sin I + r 30 


(42) 


i 3 sin I 3 

V = n cos I 


Z<f> 


3r 


r 30 r sin 6 30 


(43) 


h 1 


r kn 
I efTr 


kn 3(T e , N) 
efTr 3( r , F) ds 


i: 


3 ( N ’ T -> , 

eNr sin 0 3(s, 0) ds 


(44) 


Here we denote the Jacobian of any two functions f , g of variables C, 7 ) by 

3( f , g) _ 3f 3g 3f 3g 
V) 3 C 3 V 

The drift arising from the first integral on the right-hand side will be in the EW 
direction, while the drift arising from the second term will be in a meridian 
plane. It follows that as both N and T e are likely to vary most rapidly in the 
vertical direction, the major contribution to the continuity equation is likely to 
arise from the second term on the right-hand side of Equation (44). The order 
of magnitude of the important drift component (judged by this criterion) is 
therefore, 


kT e 

~ 3 x 10 2 cm/sec , (45) 

where T e ~ 1000°K, using Equations (3) and (41) with r ~ 7 x 10 8 cm. This also 
is negligibly small. 
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A Magnetohydrodynamic (orE x B ) Drift 

The electric field which gives rise to this drift is such that according to 
Equation (18), is constant along a field line. This particular form of electro- 
dynamic drift is well known. Martyn 8 first suggested the importance of this type 
of transport in the F2 region, and according to calculations by Moffett and Hanson 9 
and Bramley and Peart 10 it plays a great part in the geomagnetic control of the 
F2 layer. 

THE CONTINUITY EQUATION 

Denote by D a the coefficient of ambipolar diffusion defined by 

2kT i Da' 1 = \ K v in + 2m, 5 v en ) (46) 


Equation (10) then gives 




II = 


D„ 




2NkT. 


7 fe" (p e + Pi ) - Nin i e"_ 


+ 


Z*e V en 


eN 


(m. v. 

\ i in 


+ 2m 



(47) 


At this stage it is possible to show finally that may be omitted. Equations (6) 
and the values m e ~ 10 -27 gm, ~ 10 -23 gm give 


1T1 - 1 

® 12. ~ 4 x in 3 

m. v. H x AU 

i in 


(48) 


Also Dougherty’s estimate of 10 ' 12 e.m.u. for J 11 , together with II ~ 10 6 cm -3 and 
Equation (3) give 



~ 10 2 cm/sec . 


(49) 


As | v. 11 | is likely to be of order of magnitude lm/sec it appears that the last term 
in Equation (47) is likely to be inappreciable; in comparison so that for most 
purposes 


Nv," = 



(50) 
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It follows that if Q denotes the rate of production, and L the loss rate, a good ap- 
proximation to the F2 region continuity equation is 


<3N 

3t 


Q - 


L - div(Nv n 11 ) - div 



+ 


l D * 

dlv l 2kT7 


V 1 ' 


(P e + Pi) ' N mi 



(51) 


where v 11 denotes the known component of the neutral air velocity along a field 
line, and v. denotes the electrodynamic drift. 



B 

& * 8 rad Q 3 


(52) 


at right angles to a field line . 


CONCLUSION 


Under normal F2 region conditions the standard view of the ambipolar dif- 
fusion problem appears to be correct. The ion-electron pairs may be thought of 
as sliding down the magnetic field lines with a diffusive motion, while at the same 
time the whole field line drifts at right angles to itself. The continuity equation 
is then Equation (51). 

The current between conjugate points of the dynamo layer (Dougherty 1 ) 
appears to have little effect if its magnitude is as small as 10 " 7 amps/m 2 . How- 
ever, if were to become 100 times larger, its effects on v i 11 might become 
appreciable. Even then, its effects upon V • Nv, 11 would be noticeable only when 
v en / v i n in Equation (47) varies appreciably in a vertical distance of one scale 
height. (This would occur, for example, in a thin layer where the electron tem- 
perature were greater than the ion temperature.) 
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